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ABSTRACT 

Context. Large Scale Structures (LSS) in the universe can be traced using the neutral atomic hydrogen Hi through its 21 cm emission. 
Such a 3D matter distribution map can be used to test the Cosmological model and to constrain the Dark Energy properties or its 
equation of state. A novel approach, called intensity mapping can be used to map the Hi distribution, using radio interferometers with 
large instanteneous field of view and waveband. 

Aims. In this paper, we study the sensitivity of different radio interferometer configurations, or multi-beam instruments for the 
observation of large scale structures and BAO oscillations in 21 cm and we discuss the problem of foreground removal. 
Methods. For each configuration, we determine instrument response by computing the (u,v) or Fourier angular frequency plane 
coverage using visibilities. The (u,v) plane response is the noise power spectrum, hence the instrument sensitivity for LSS P(k;) 
measurement. We describe also a simple foreground subtraction method to separate LSS 21 cm signal from the foreground due to the 
galactic synchrotron and radio sources emission. 

Results. We have computed the noise power spectrum for different instrument configuration as well as the extracted LSS power 
spectrum, after separation of 21cm-LSS signal from the foregrounds. We have also obtained the uncertainties on the Dark Energy 
parameters for an optimized 21 cm BAO survey. 

Conclusions. We show that a radio instrument with few hundred simultaneous beams and a collecting area of ~ lOOOOm^ will be 
able to detect BAO signal at redshift z ~ 1 and will be competitive with optical surveys. 

Key words, large-scale structure of Universe - dark energy - Instrumentation: interferometers - Radio lines; galaxies - Radio 
continuum: general 



1. Introduction 

The study of the statistical properties of Large Scale Structure 
(LSS) in the Universe and their evolution with redshift is 
one the major tools in observational cosmology. These struc- 
tures are usually mapped through optical observation of galax- 
ies which are used as a tracer of the underlying matter 
distribution. An alternative and elegant approach for map- 
ping the matter distribution, using neutral atomic hydrogen 
(Hi ) as a tracer with intensity mappin g , has been proposed 
in recent years ( [Peterson et al. (2006)| ) ( [Chang et al. (2008) i. 
Mapping the matter distribution using HI 21 cm emis- 
sion as a tracer has been extensively discussed in literature 
( [Furlanetto et al. (2006)| ) ( [Tegmark & Zaldarriaga (2"009)]l and is 
being used in projects such as LOFAR ('Rotter ing et al. (2006)| l 
or MWA (Bowman et al. (2009) ) to observe reionisation at red- 
shifts z ~ 10. 

Evidence in favor of the acceleration of the expansion 
of the universe have been accumulated over the last twelve 
years, thanks to the observation of distant supernovae, CMB 
anisotropics and detailed analysis of the LSS. A cosmological 
Constant (A) or new cosmological energy density called Dark 
Energy has been advocated as the origin of this acceleration. 



Dark Energy is considered as one of the most intriguing puzzles 
in Physics and Cosmology. Several cosmological probes can be 
used to constrain the properties of this new cosmic fluid, more 
precisely its equation of state: The Hubble Diagram, or luminos- 
ity distance as a function of redshift of supernovae as standard 
candles, galaxy clusters, weak shear observations and Baryon 
Acoustic Oscillations (BAO). 

BAO are features imprinted in the distribution of galax- 
ies, due to the frozen sound waves which were present 
in the photon-baryon plasma prior to recombination at z 
~ 1100. This scale can be considered as a standard ruler 
with a comoving length of ~ ISOMpc. These features have 
been first observed in the CMB anisotropies and are usu- 
ally referred to as acoustic peaks ( [Mauskopf et al. (2000) 



[Larson et al. (201 l)l l. The BAO modulation has been subse 
quently observed in the distribution of galaxies at low red- 
shift { z < 1) in the galaxy-galaxy correlation function 
by the SDSS ([Eisenstein et al. (2005)]l ([Percival et al. (2007)|l 



( [Percival et al. (2010)|l, 2dGF RS ( |Cole et al. (2005)[ ) as weU as 
WiggleZ ( Blake et al. (2011)l l optical galaxy surveys. 

Ongoing ( [Eisenstein et al. (201 1)[ ) or future surveys 
(ILSST.Sciencel) plan to measure precisely the BAO 
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scale in the redshift range < z < 3, using either 
optical observation of galaxies or through 3D map- 
ping Lyman a absorption lines toward distant quasars 
( [McDonald et al. (2006)^ , ( [McDonald & Eisenstein (2007)l l. 
Radio observation of the 21 cm emission of neutral hydrogen 
appears as a very promising technique to map matter distribu- 
tion up to redshift z ~ 3, complementary to optical surveys, 
especially in the optical redshift desert range 1 < z < 2, and 



possibly up to the reionization redshift ( Wyithe et al.(2008) i. 

In section 2, we discuss the intensity mapping and its poten- 
tial for measurement of the Hi mass distribution power spectrum. 
The method used in this paper to characterize a radio instrument 
response and sensitivity for P\i^{k) is presented in section 3. We 
show also the results for the 3D noise power spectrum for sev- 
eral instrument configurations. The contribution of foreground 
emissions due to the galactic synchrotron and radio sources is 
described in section 4, as well as a simple component separation 
method. The performance of this method using two different sky 
models is also presented in section 4. The constraints which can 
be obtained on the Dark Energy parameters and DETF figure of 
merit for typical 21 cm intensity mapping survey are discussed 
in section 5. 

2. Intensity mapping and Hi power spectrum 

2.1. 21 cm intensity mapping 

Most of the cosmological information in the LSS is located at 
large scales (> Ideg), while the interpretation at smallest scales 
might suffer from the uncertainties on the non linear clustering 
effects. The BAO features in particular are at the degree angular 
scale on the sky and thus can be resolved easily with a rather 
modest size radio instrument (diameter D < 100 m). The spe- 
cific BAO clustering scale (^bao) can be measured both in the 
transverse plane (angular correlation function, (k-^^o) or along 

the longitudinal (line of sight or redshift (k^^^Q) direction. A di- 
rect measurement of the Hubble parameter H{z) can be obtained 
by comparing the longitudinal and transverse BAO scales. A rea- 
sonably good redshift resolution 6z < 0.01 is needed to resolve 
longitudinal BAO clustering, which is a challenge for photomet- 
ric optical surveys. 

In order to obtain a measurement of the LSS power spectrum 
with small enough statistical uncertainties (sample or cosmic 
variance), a large volume of the universe should be observed, 
typically few Gpc"*. Moreover, stringent constraint on DE pa- 
rameters can only be obtained when comparing the distance or 
Hubble parameter measurements with DE models as a function 
of redshift, which requires a significant survey depth Az > 1 . 

Radio instruments intended for BAO surveys must thus have 
large instantaneous field of view (FOV > lOdeg^) and large 
bandwidth (Av > 100 MHz) to explore large redshift domains. 

Although the application of 21 cm radio survey to cosmol- 
ogy, in particular LSS mapping has been discussed in length in 
the framework of large future instruments, such as the SKA (e.g 
ISKA.Scien"cel [Abdalla & Rawhngs (2005)[ ), the method envis- 
aged has been mostly through the detection of galaxies as Hi 
compact sources. However, extremely large radio telescopes are 
required to detected Hi sources at cosmological distances. The 
sensitivity (or detection threshold) limit S for the total power 
from the two polarisations of a radio instrument characterized by 
an effective collecting area A, and system temperature Ts^s can 
be written as 
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Table 1. Sensitivity or source detection limit for 1 day integra- 
tion time (86400 s) and 1 MHz frequency band (left). Source 21 
cm brightness for IO'^Mq Hi for different redshifts (right) 



where f„„ is the total integration time and Sv is the detection fre- 
quency band. In table [T] (left) we have computed the sensitivity 
for 6 different sets of instrument effective area and system tem- 
perature, with a total integration time of 86400 seconds (1 day) 
over a frequency band of 1 MHz. The width of this frequency 
band is well adapted to detection of Hi source with an intrin- 
sic velocity dispersion of few 100 km/s. These detection limits 
should be compared with the expected 21 cm brightness 521 of 
compact sources which can be computed using the expression 
below (e.g |Binney & Merrifield (1998)| l : 



>21 



Mh, 

0.021/iJy — X 



i^)x?«^(l.z) (2) 



where Mh, is the neutral hydrogen mass, c/l(z) is the luminosity 
distance and cr,, is the source velocity dispersion. 

In table[T](right), we show the 21 cm brightness for compact 
objects with a total Hi mass of IO'^Mq and an intrinsic veloc- 
ity dispersion of 200 km/s. The luminosity distance is computed 
for the standard WMAP ACDM universe. 10*^ - 10'°Mq of n eu- 
tral gas mass is typical for large galaxies ( [Lah et al. (2009)| l. It 
is clear that detection of Hi sources at cosmological distances 
would require collecting area in the range of lO^m^. 

Intensity mapping has been suggested as an alternative 
and economic method to map the 3D distribution of neu- 
tral hydrogen by ([Chang et al. (2008) i and further studied by 
( [Ansari et al. (2008)| l ( [Seo et al (20T0)| l. In this approach, sky 
brightness map with angular resolution ~ 10-30 arc.min is made 
for a wide range of frequencies. Each 3D pixel (2 angles 0, fre- 
quency V or wavelength A) would correspond to a cell with a 
volume of ~ lO^Mpc-*, containing ten to hundred galaxies and a 
total Hi mass ~ IO'^Mq. If we neglect local velocities relative to 
the Hubble flow, the observed frequency v would be translated 
to the emission redshift z through the well known relation: 

z(v) = Y^lZl ■ v(2) = with V2I = 1420.4 MHz (3) 

(l+z) 

; A{z) = /I21 X (1 -H z) with A21 = 0.21 1 m (4) 
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The large scale distribution of the neutral hydrogen, down to an- 
gular scales of ~ lOarc.min can then be observed without the 
detection of individual compact Hi sources, using the set of sky 
brightness map as a function of frequency (3D-brightness map) 
B2i{&, A). The sky brightness B21 (radiation power/unit solid an- 
gle/unit surface/unit frequency) can be converted to brightness 
temperature using the well known black body Rayleigh-Jeans 
approximation: 

B(T,A) = ^ 

A- 



2 



R. Ansari et al.: 21 cm observation of LSS at z ~ 1 



2.2. Hi power spectrum and BAO 

In the absence of any foreground or background radiation, the 
brightness temperature for a given direction and wavelength 
T2i(&, A) would be proportional to the local Hi number density 
«H/(€), z) through the relation: 



3 h 
327r kg 



All X 



H(z) 



{l+zfxnHi(&,z) (5) 



where A21 = 2.85 10^'^s"' ( |Lang (1999)) is the spontaneous 21 
cm emission coefficient, h is the Planck constant, c the speed of 
light, kg the Boltzmann constant and H{z) is the Hubble parame- 
ter at the emission redshift. For a ACDM universe and neglecting 
radiation energy density, the Hubble parameter can be expressed 
as: 

H(z) ^ hm [a„(l + zf + Qa]' X 100 km/s/Mpc (6) 

Introducing the Hi mass fraction relative to the total baryon mass 
fn,, the neutral hydrogen number density relative fluctuations 
can be written as, and the corresponding 21 cm emission tem- 
perature can be written as: 



/h,(z)xQ,^x(^(0,z) + 1 
mn \ Ph, 

6pH. 



(7) 



(8) 



r2i(0,^(z)) = r2i(z)x^(0,z) + l 
\Ph, 

where Q.B,Pcm are respectively the present day mean baryon cos- 
mological and critical densities, nin is the hydrogen atom mass, 

and is the Hi density fluctuations. 

The present day neutral hydrogen fraction /h,(0) present in 
local galaxies has been measured to be ~ 1% of the baryon den- 
sity ( |Zwaan et aL(20 05)): 



Q.H, ^ 3.5 10" 



0.008 X Qb 



The neutral hydrogen fraction is expected to increase with red- 
shift, as gas is used in star formation during galaxy formation 
and evolution. Study of Lyman-a absorption indicate a factor 3 
increase in the neutral hydrogen fraction at z = 1.5 in the in- 
tergalactic medium ( |Wolf et al.(2 005) ), compared to its present 
day value /h,(z - 1.5) ~ 0.025. The 21 cm brightness temper- 
ature and the corresponding power spectrum can be written as 
( IBarkana & Loeb (2007)| and |Madau et al. (1997)| ) : 



PT.Ak) = [flliz)) P(k) 

T2i(z) ^ 0.084 mK 



(1 +zfhm) 
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(9) 
(10) 



The table |2] shows the mean 21 cm brightness tempera- 
ture for the standard ACDM cosmology and either a constant 
Hi mass fraction fn, - 0.01, or linearly increasing fn, - 
0.008 X (1 + z). Figure[T]shows the 21 cm emission power spec- 
trum at several redshifts, with a constant neutral fraction at 2% 
(/h; = 0.02). The matter power spectrum has been computed us- 
ing the |Eisenstein & Hu (1998)] parametrisation. The correspon- 
dence with the angular scales is also shown for the standard 
WMAP ACDM cosmology, according to the relation: 
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where k is the comoving wave vector and c/a(z) is the angular 
diameter distance. 
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Table 2. Mean 21 cm brightness temperature in mK, as a func- 
tion of redshift, for the standard ACDM cosmology with con- 
stant Hi mass fraction at //f,(z)=0.01 (a) or linearly increasing 
mass fraction (b) fyXz) = 0.008(1 + z) 
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Fig. 1. Hi 21 cm emission power spectrum at redshifts z-\ (blue) 
and z=2 (red), with neutral gas fraction fn, - 2% 

3. interferometric observations and P(k) 
measurement sensitivity 

3.1. Instrument response 



We introduce briefly here the principles of interferomet 
ric observations and the definition of quantities 

may 



ful for our calculations. Interested reader 



use- 
refer to 

(Thompson, Moran & S wenson (2001)| l for a detailed and com- 
plete presentation of observation methods and signal processing 
in radio astronomy. In astronomy we are usually interested in 
measuring the sky emission intensity, I{&,A) in a given wave 
band, as a function of the sky direction. In radio astronomy 
and interferometry in particular, receivers are sensitive to the 
sky emission complex amplitudes. However, for most sources, 
the phases vary randomly with a spatial correlation length 
significantly smaller than the instrument resolution. 

I(Q,A) = \Ai&,A)\^ , I eR,AeC (12) 
< A{&,A)A*{0',A) >,ime^ for 5t 0' (13) 

A single receiver can be characterized by its angular complex 
amplitude response B(0, v) and its position r in a reference 
frame, the waveform complex amplitude s measured by the re- 
ceiver, for each frequency can be written as a function of the 
electromagnetic wave vector ^£jvf(0. A) : 



s(A) 



If 



dQ A(0, A)B(Q, A)e' 



i(kEM-r) 



(14) 



We have set the electromagnetic (EM) phase origin at the center 
of the coordinate frame and the EM wave vector is related to the 
wavelength A through the usual equation I^^mI = 2n/A. The re- 
ceiver beam or antenna lobe L(0, A) corresponds to the receiver 
intensity response: 

L(0, A) = B(0, A) B*{&, A) (15) 
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The visibility signal of two receivers corresponds to the time av- 
eraged correlation between signals from two receivers. If we as- 
sume a sky signal with random uncorrected phase, the visibility 
"Vn signal from two identical receivers, located at the position 
ri and can simply be written as a function of their position 
difference \r - ri - 



'VniA) si{A)s2{Ar > 



-ff 



dQ> 7(0, A)L(&, ^)e'(*™ '^'-) 



(16) 

This expression can be simplified if we consider receivers with 
narrow field of view (L{0,A) ^ for |0| > lOdeg. ), and 
coplanar in respect to their common axis. If we introduce two 
Cartesian like angular coordinates {a,/3) centered at the com- 
mon receivers axis, the visibilty would be written as the 2D 
Fourier transform of the product of the sky intensity and the re- 
ceiver beam, for the angular frequency (m, v)i2 = 27r( 



"VniA) 



ff 



dad/3 I(a,f3)L(a,/3)exp 



i2n 



Ax Ay 



A 



A 



(17) 

where (Ax, Ay) are the two receiver distances on a plane perpen- 
dicular to the receiver axis. The x and y axis in the receiver plane 
are taken parallel to the two ia,/3) angular planes. 

Furthermore, we introduce the conjugate Fourier variables 
(m, v) and the Fourier transforms of the sky intensity and the re- 
ceiver beam: 



I(a,/3, A) 
L{a,l3, A) 



(u,v) 
I(u, V, A) 
XXu, V, A) 



The visibility can then be interpreted as the weighted sum 
of the sky intensity, in an angular wave number domain located 
around (m, v)i2 = 27r(^, ^). The weight function is given by the 
receiver beam Fourier transform. 

^n(A)- JJ dudv I{u,v, A) Ziu- 271^, v-2n^. A) (18) 



A 



A single receiver instrument would measure the total power 
integrated in a spot centered around the origin in the (m, v) or the 
angular wave mode plane. The shape of the spot depends on the 
receiver beam pattern, but its extent would be ~ 27tD/A, where 
D is the receiver physical size. 

The correlation signal from a pair of receivers would mea- 
sure the integrated signal on a similar spot, located around the 
central angular wave mode (m, v)i2 determined by the relative 
position of the two receivers (see figure |2|i. In an interferometer 
with multiple receivers, the area covered by different receiver 
pairs in the (m, v) plane might overlap and some pairs might 
measure the same area (same base lines). Several beams can be 
formed using different combination of the correlations from a set 
of antenna pairs. 

An instrument can thus be characterized by its (m, v) plane 
coverage or response Kiu, v. A). For a single dish with a single 
receiver in the focal plane, the instrument response is simply the 
Fourier transform of the beam. For a single dish with multiple 
receivers, either as a Focal Plane Array (FPA) or a multi-horn 
system, each beam (b) will have its own response 'Rbiu, v. A). 
For an interferometer, we can compute a raw instrument re- 
sponse 'Rrawiu, V, A) which corresponds to (u, v) plane coverage 
by all receiver pairs with uniform weighting. Obviously, differ- 
ent weighting schemes can be used, changing the effective beam 
shape and thus the response '??,v(m, v. A) and the noise behaviour 
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Fig. 2. Schematic view of the (u, v) plane coverage by interfero- 
metric measurement. 



If the same Fourier angular frequency mode is measured by sev- 
eral receiver pairs, the raw instrument response might then be 
larger that unity. This non normalized instrument response is 
used to compute the projected noise power spectrum in the fol- 



lowing section (3.3 i. We can also define a normalized instrument 
response, 'R„o„„(m, v,A)< 1 as: 

KonJu, V, A) = :R(m, V, ^)/MaX(u,v) YRiu, V, A)] (19) 

This normalized instrument response can be used to compute the 
effective instrument beam, in particular in section |43 



3.2. Noise power spectrum 

Let's consider a total power measurement using a receiver at 
wavelength A, over a frequency bandwidth 6v centered on vq, 
with an integration time f,,,,, characterized by a system temper- 
ature Tsys. The uncertainty or fluctuations of this measurement 

due to the receiver noise can be written as cr^ . = ^-7^. This 

noise tin, ov 

term corresponds also to the noise for the visibility 'Vn mea- 
sured from two identical receivers, with uncorrected noise. If 
the receiver has an effective area A ^ ttD^/A- or A ^ D^Dy, 
the measurement corresponds to the integration of power over a 
spot in the angular frequency plane with an area ~ AjA^. The 
noise spectral density, in the angular frequencies plane (per unit 
area of angular frequencies || x |j), corresponding to a visibility 
measurement from a pair of receivers can be written as: 



rjpair 



tin,5v D2 



units : x rad^ 



(20) 
(21) 



The sky temperature measurement can thus be characterized 
by the noise spectral power density in the angular frequencies 



plane P 



in Kelvin per unit area of angular frequen- 



2k ^ lir' '^'^^ characterize the sky temperature measure- 
ment with a radio instrument by the noise spectral power density 
in the angular frequencies plane Pnoiseiu, v) in units of Kelvin^ 
per unit area of angular frequencies x For an interfer- 
ometer made of identical receiver elements, several (n) receiver 
pairs might have the same baseline. The noise power density in 
the corresponding (u, v) plane area is then reduced by a factor 
More generally, we can write the instrument noise spectral 
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power density using the instrument response defined in section 
[TTI : 

PnoisAu, V) = 7^^^^^ (22) 

When the intensity maps are projected in a three dimensional 
box in the universe and the 3D power spectrum P(k) is com- 
puted, angles are translated into comoving transverse distances, 
and frequencies or wavelengths into comoving radial distance, 
using the following relations: 



6v 

6u,6v 
1 

Sv 



c 6v 



6kn - 



H{z) V 
6u , 6v 

(1 +z)dA(z) 
H(z) 1 V 



(1+z) 



H{z) 



-5v 



H(z) 1 V21 



(1 + z) dv c (1 +z)^ 6v 



plane corresponding to the equation 2 1 above, the three dimen 



sional projected noise spectral density can then be written as: 



P noise^K) ~ 2 



sys 



dliz) 



Hiz) 



(l+z)' 



(27) 



Piwiseih) would be in units of mK^ Mpc-' with Tsys expressed 
in mK, f,„, is the integration time expressed in second, V21 in Hz, 
c in km/s, c/^ in Mpc and H{z) in km/s/Mpc. 

The matter or Hj distribution power spectrum determination 
statistical errors vary as the number of observed Fourier modes, 
which is inversely proportional to volume of the universe which 
is observed (sample variance). As the observed volume is pro- 
portional to the surveyed solid angle, we consider the survey of 
a fixed fraction of the sky, defined by total solid angle Q.,ot, per- 
formed during a determined total observation time tgbs- A single 
dish instrument with diameter D would have an instantaneous 
field of view 0.fov ~ {j^ ^ ™d would require a number of point- 
ings N point — to cover the survey area. Each sky direction 
or pixel of size tlfov will be observed during an integration time 
tint - tohsIN point- Using equation 27 and the previous expression 



for the integration time, we can compute a simple expression for 
the noise spectral power density by a single dish instrument of 
diameter D: 



r^survey 



(k) = 2 



tobs V2I 



dA(z)j^(l+zr (28) 



It is important to note that any real instrument do not have 
a flat response in the (m, v) plane, and the observations pro- 
vide no information above a certain maximum angular frequency 
Umax, Vmax- One has to take into account either a damping of the 
observed sky power spectrum or an increase of the noise spectral 
power if the observed power spectrum is corrected for damping. 
The white noise expressions given below should thus be con- 
sidered as a lower limit or floor of the instrument noise spectral 
density. 

For a single dish instrument of diameter D equipped with a 
multi-feed or phase array receiver system, with independent 
beams on sky, the noise spectral density decreases by a factor N, 
thanks to the increase of per pointing integration time: 



2 Tls^t, 



■ dliz) 



Hiz) 



(29) 
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Fig. 3. Minimal noise level for a 100 beams instrument with 
Tsys - 50K as a function of redshift (top). Maximum k value for 
a 100 meter diameter primary antenna (bottom) 



This expression (eq. 29 1 can also be used for a filled interfer- 
ometric array of N identical receivers with a total collection area 
~ D^. Such an array could be made for example of N = q X q 
small dishes, each with diameter D/q, arranged as, qx q square. 

For a single dish of diameter D, or an interferometric instru- 
ment with maximal extent D, observations provide information 



up to u„ 



2nD/A. This value of u„ 



would be 



mapped to a maximum transverse cosmological wave number 
k^ ■ 

max' 

^"'^^ kL.<-^^ (30) 



k^ = 



(l+z)dA dA(l+z)^A2i 

Figure [3] shows the evolution of the noise spectral density 
f """"VA:) as a function of redshift, for a radio survey of the sky, 

noise ^ ^ ' J J ■> 

using an instrument with A = 100 beams and a system noise 
temperature r,, ^ - 50K. The survey is supposed to cover a quar- 
ter of sky Q.tot - srad, in one year. The maximum comoving 
wave number kmax is also shown as a function of redshift, for an 
instrument with D = 100 m maximum extent. In order to take 
into account the radial component of k and the increase of the 
instrument noise level with fc^, we have taken the effective k,„ax 



as half of the maximum transverse k^^^^. of eq. 30 

^maxiz) — 



D = 100m 



dAd+z)'^ A21 



(31) 



3.3. Instrument configurations and noise power spectrum 

We have numerically computed the instrument response 
Kiu, V, A) with uniform weights in the (u, v) plane for several in- 
strument configurations: 

a : A packed array of n = 121 Dj/s/, = 5 m dishes, arranged in a 
square 11x11 configuration (q = 11). This array covers an 
area of 55 x 55 m^ 
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(b) 128 D=5m dishes in 8 rows 



(c) 129 D=5m dishes 
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Fig. 4. Array layout for configurations (b) and (c) with 128 and 
129 D=5 meter diameter dishes. 



b : An array of n - 128 Dj/j/, = 5 m dishes, arranged in 8 rows, 
each with 16 dishes. These 128 dishes are spread over an area 
80 X 80 m^. The array layout for this configuration is shown 
in figure |4] 

c : An array of « = 129 Ddish = 5 m dishes, arranged over an 
area 80 x 80 m^. This configuration has in particular 4 sub- 
arrays of packed 16 dishes (4 x 4), located in the four array 
corners. This array layout is also shown figure|4] 

d : A single dish instrument, with diameter D = 75 m, equipped 
with a 100 beam focal plane receiver array. 

e : A packed array of n = 400 Ddish = 5 m dishes, arranged in a 
square 20 x 20 configuration (q - 20). This array covers an 
area of 100 x 100 m^ 

f : A packed array of 4 cylindrical reflectors, each 85 meter 
long and 12 meter wide. The focal line of each cylinder is 
equipped with 100 receivers, each 2A long, corresponding to 
~ 0.85 m at z = 1. This array covers an area of 48 x 85 m^, 
and have a total of 400 receivers per polarisation, as in the (e) 
configuration. We have computed the noise power spectrum 
for perfect cylinders, where all receiver pair correlations are 
used (fp), or for a non perfect instrument, where only corre- 
lations between receivers from different cylinders are used. 

g : A packed array of 8 cylindrical reflectors, each 102 meter 
long and 12 meter wide. The focal line of each cylinder is 
equipped with 120 receivers, each 2A long, corresponding to 
~ 0.85 m at z = 1. This array covers an area of 96 x 102 m^ 
and has a total of 960 receivers per polarisation. As for the (f) 
configuration, we have computed the noise power spectrum 
for perfect cylinders, where all receiver pair correlations are 
used (gp), or for a non perfect instrument, where only corre- 
lations between receivers from different cylinders are used. 

We have used simple triangular shaped dish response in the 
(m, v) plane. However, we have introduced a filling factor or il- 
lumination efficiency 77, relating the effective dish diameter Dm 
to the mechanical dish size D'" - rj Ddish- The effective area 
Ae oc rp- scales as rj^ or rj^riy. 



[±1jiDII'IA] 



sin(;r(D'"/^) sin0) 



7r{D'"/A) sine 



For the multi-dish configuration studied here, we have taken the 
illumination efficiency factor 77 = 0.9. 

For the receivers along the focal line of cylinders, we have 
assumed that the individual receiver response in the (m, v) plane 
corresponds to one from a rectangular shaped antenna. The il- 
lumination efficiency factor has been taken equal to 77^ - 0.9 



in the direction of the cylinder width, and rjy - 0.8 along the 
cylinder length. It should be noted that the small angle approxi- 
mation used here for the expression of visibilities is not valid for 
the receivers along the cylinder axis. However, some preliminary 
numerical checks indicate that the results obtained here for the 
noise spectral power density would not change significantly. The 
instrument responses shown here correspond to fixed pointing 
toward the zenith, which is the case for a transit type telescope. 



^□(M,V,^)= l\ {U)X l\ (V) 

[±2nD'i'lA] [±2nD'l}lA] 



(34) 



Figure [5] shows the instrument response Hiu, v. A) for the four 
configurations (a,b,c,d) with ~ 100 receivers per polarisation. 
The resulting projected noise spectral power density is shown in 
figure |6] The increase of P„oiAk) at low k"'"""' < 0.02 is due to 
the fact that we have ignored all auto-correlation measurements. 
It can be seen that an instrument with 100-200 beams and Tsys - 
50K should have enough sensitivity to map LSS in 21 cm at 
redshift z=l. 



4. Foregrounds and Component separation 

Reaching the required sensitivities is not the only difficulty of 
observing the large scale structures in 21 cm. Indeed, the syn- 
chrotron emission of the Milky Way and the extra galactic ra- 
dio sources are a thousand times brighter than the emission 
of the neutral hydrogen distributed in the universe. Extracting 
the LSS signal using Intensity Mapping, without identifying 
the Hi point sources is the main challenge for this novel ob- 
servation method. Although this task might seem impossible 
at first, it has been suggested that the smooth frequency de- 
pendence of the synchrotron emissions can be used to separate 
the faint LSS signal from the Galactic and radio source emis- 
sions. However, any real radio instrument has a beam shape 
which changes with frequency: this instrumental effect signifi- 
cantly increases the difficulty and complexity of this component 
separation technique. The effect of frequency dependent beam 
shape is some time referred to as mode mixing. See for example 
(Morales et al. (2006)| l, (Bowman et al. (2009)| l. 

In this section, we present a short description of the fore- 
ground emissions and the simple models we have used for com- 
puting the sky radio emissions in the GHz frequency range. We 
present also a simple component separation method to extract 
the LSS signal and its performance. We show in particular the 
effect of the instrument response on the recovered power spec- 
trum. The results presented in this section concern the total sky 
emission and the LSS 21 cm signal extraction in the z ~ 0.6 
redshift range, corresponding to the central frequency v ~ 884 
MHz. 



(32) 4.1. Synchrotron and radio sources 

(33) 



We have modeled the radio sky in the form of three dimensional 
maps (data cubes) of sky temperature brightness T{a, 6, v) as a 
function of two equatorial angular coordinates (a, S) and the fre- 
quency V. Unless otherwise specified, the results presented here 
are based on simulations of 90 x 30 - 2500 deg^ of the sky, 
centered on a = lOhOOm, (5 = +10deg., and covering 128 MHz 
in frequency. We have selected this particular area of the sky in 
order to minimize the Galactic synchrotron foreground. The sky 
cube characteristics (coordinate range, size, resolution) used in 
the simulations are given in the table |3] 
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u-v coverage, (a) 11 x11 D=5m dishes Over 55mx55m 



u-v coverage, (b) 8 row of 1 6xD=5m dishes Over SOmxSOm 



'PIOOO 



500 



-500 - 



-1000 




-1000 -500 500 1000 

u (Radian'-I) u=1000->~21 arcmin 

u-v coverage, (c) 1 29 D=5m dishes Over SOmxSOm 



91000 



500 



-500 - 



-1000 




'PIOOO 



500 



-500 



-1000 




-1000 -500 500 1000 

u (Radian'^-1) u=1000->~21 arcmin 

u-v coverage, (d) D=75 m dish with 1 00 beams 




-1000 -500 500 1000 

u (Radian'-I) u=1000->~21 arcmin 



-1000 -500 500 1000 

u (Radian'-1) u=1000->~21 arcmin 



Fig. 5. (u,v) plane coverage (raw instrument response 'R(u,v,A) for four configurations, (a) 121 Z)^,v, = 5 meter diameter dishes 
arranged in a compact, square array of 1 1 x 1 1, (b) 128 dishes arranged in 8 row of 16 dishes each (fig.|4]l, (c) 129 dishes arranged 
as shown in figure]?], (d) single D=75 meter diameter, with 100 beams, (color scale : black < 1, blue, green, yellow, red > 80) 



Two different methods have been used to compute the sky 
temperature data cubes. We have used the Global Sky Model 
(GSM) ( [Ohveira-Costa et al. (2008)| l tools to generate full sky 
maps of the emission temperature at diff'erent frequencies, from 
which we have extracted the brightness temperature cube for the 
region defined above (Model-I/GSM Tg^mia, 6, v)). As the GSM 
maps have an intrinsic resolution of ~ 0.5 degree, it is difficult 
to have reliable results for the effect of point sources on the re- 
constructed LSS power spectrum. 

We have thus made also a simple sky model us- 
ing the Haslam Galactic synchrotron map at 408 MHz 
dHaslam et al. (1982)D and the N RAO VLA Sky Surve y (NVSS) 
1.4 GHz radio source catalog ( Condon et al. (1998)| l. The sky 
temperature cube in this model (Model-II/HaslamH-NVSS) has 
been computed through the following steps: 



1. The Galactic synchrotron emission is modeled as a power 
law with spatially varying spectral index. We assign a power 
lawindex/3 - -2.8+0.15 to each sky direction. y6 has a gaus- 
sian distribution centered at -2.8 and with standard deviation 
cTg = 0.15. The synchrotron contribution to the sky tempera- 
ture for each cell is then obtained through the formula: 

2. A two dimensional T^yssio, 6) sky brightness temperature at 
1.4 GHz is computed by projecting the radio sources in the 
NVSS catalog to a grid with the same angular resolution as 
the sky cubes. The source brightness in Jansky is converted 
to temperature taking the pixel angular size into account (~ 
21mK/mJansky at 1.4 GHz and 3' x 3' pixels). A spectral 
index fSsrc £ [-1-5, -2] is also assigned to each sky direction 



7 



R. Ansari et al.: 21 cm observation of LSS at z ~ 1 



PNoise(k) : 121,128,129,100 receiver/polar @z=1 



PNoise(k) : Dishes/Cylinders, 400/400/960 recv/pol @z=1 




2000 



0.006 0.01 0.02 0.04 0.06 0.1 0.2 0.4 

k_comov (h Mpc'^-I) 




0.006 0.01 0.02 0.04 0.06 0.1 0.2 0.4 

k_comov (h Mpc^-I) 



Fig. 6. P(k) LSS power and noise power spectrum for several interferometer configurations ((a),(b),(c),(d),(e)Xf)Xg)) with 121, 128, 
129, 400 and 960 receivers. 





range 


center 


Right ascension 


105 < a < 195 deg. 


150 deg. 


Declination 


-5 < 5 < 25 deg. 


+ 10 deg. 


Frequency 


820 < V < 948 MHz 


884 MHz 


Wavelength 


36.6 < A < 31.6 cm 


33.9 cm 


Redshift 


0.73 < z < 0.5 


0.61 




resolution 


N-cells 


Right ascension 


3 arcmin 


1800 


Declination 


3 arcmin 


600 


Frequency 


500 kHz {dz -10-3) 


256 



Table 3. Sky cube characteristics for the simulation performed 
in this paper Cube size : 90 deg. x 30 deg. x 128 MHz 1800 x 
600 X 256 123 10^ cells 



for the radio source map; we have taken /3src as a flat random 
number in the range [-1.5,-2], and the contribution of the 
radiosources to the sky temperature is computed as follows: 



^rat/ire (Q', ^, v) — Tjjyss X 



1 1420 MHz/ 



3. The sky brightness temperature data cube is obtained 
through the sum of the two contributions. Galactic syn- 
chrotron and resolved radio sources: 

Tfg„d(a, 6, v) = T,y„c(a, 6, v) + Tradsrda, 6, v) 

The 21 cm temperature fluctuations due to neutral hydrogen 
in large scale structures Tiss{a, S, v) have been computed using 
the SimLSS[^software package: complex normal Gaussian fields 
were first generated in Fourier space. The amplitude of each 



mode was then multiplied by the square root of the power spec- 
trum P{k) at z = computed according to the parametrization of 
( |Eisenstein & Hu (1998)| l. We have used the standard cosmolog- 
ical parameters, Hq = 71km/s/Mpc, Q,„ = 0.27, Q/, = 0.044, 
= 0.73 and w = -1. An inverse FFT was then performed 
to compute the matter density fluctuations Sp/p in the linear 
regime, in a box of 3420 x 1140 x 716 Mpc-' and evolved to 
redshift z - 0.6. The size of the box is about 2500 deg^ in the 
transverse direction and Az - 0.23 in the longitudinal direction. 
The size of the cells is 1.9 x 1.9 x 2.8 Mpc^, which correspond 
approximately to the sky cube angular and frequency resolution 
defined above. 

The mass fluctuations has been converted into temperature 
through a factor 0.13 mK, coiTesponding to a hydrogen fraction 
0.008 X (1 -H 0.6), using equation 10 The total sky brightness 
temperature is then computed as the sum of foregrounds and the 
LSS 21 cm emission: 



^sky sync ^radsrc '^Iss OR T'sky '^gsm -^i 



SimLSS : http : //www . sophya . org/SimLSS 



Iss 



Table |4] summarizes the mean and standard deviation of the 
sky brightness temperature T(a,6,v) for the different compo- 
nents computed in this study. It should be noted that the standard 
deviation depends on the map resolution and the values given in 
table |4] correspond to sky cubes computed here, with ~ 3 arc 
minute angular and 500 kHz frequency resolutions (see table|3]l. 
Figure [8] shows the comparison of the GSM temperature map at 
884 MHz with Haslam+NVSS map, smoothed with a 35 arcmin 
gaussian beam. Figure |7] shows the comparison of the sky cube 
temperature distribution for Model-I/GSM and Model-II. There 
is good agreement between the two models, although the mean 
temperature for Model-II is slightly higher (~ 10%) than Model- 
I. 
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mpfin 


stH Hpv fK'l 


Haslam 


2.17 


0.6 


NVSS 


0.13 


7.73 


Haslam+NVSS 


2.3 


7.75 


(Haslam+NVSS)*Lobe(35') 


2.3 


0.72 


GSM 


2.1 


0.8 



Table 4. Mean temperature and standard deviation for the differ- 
ent sky brightness data cubes computed for this study (see table 
[3] for sky cube resolution and size). 



BlackiGSM Red:Haslam+NVSS*Lobe(35arcmin) 




1 2 3 4 5 

SkyCube pixel temperature distribution (K) 



Fig. 7. Comparison of GSM (black) Model-II (red) sky cube 
temperature distribution. The Model-II (Haslam-i-NVSS), has 
been smoothed with a 35 arcmin gaussian beam. 

we have computed the power spectrum for the 21cm-LSS 
sky temperature cube, as well as for the radio foreground tem- 
perature cubes obtained from the two models. We have also 
computed the power spectrum on sky brightness temperature 
cubes, as measured by a perfect instrument having a 25 arcmin 
(FWHM) gaussian beam. The resulting computed power spec- 
tra are shown on figure|9] The GSM model has more large scale 
power compared to our simple Haslam-i-NVSS model, while it 
lacks power at higher spatial frequencies. The mode mixing due 
to frequency dependent response will thus be stronger in Model- 
II (Haslam+NVSS) case. It can also be seen that the radio fore- 
ground power spectrum is more than ~ 10^ times higher than the 
21 cm signal from large scale structures. This corresponds to the 
factor ~ 10^ of the sky brightness temperature fluctuations (~ 
K), compared to the mK LSS signal. 

It should also be noted that in section 3, we presented the 
different instrument configuration noise levels after correcting 
or deconvolving the instrument response. The LSS power spec- 
trum is recovered unaffected in this case, while the noise power 
spectrum increases at high k values (small scales). In practice, 
clean deconvolution is difficult to implement for real data and the 
power spectra presented in this section are NOT corrected for the 
instrumental response. The observed structures have thus a scale 
dependent damping according to the instrument response, while 
the instrument noise is flat (white noise or scale independent). 

4.2. Instrument response and LSS signal extraction 

The observed data cube is obtained from the sky brightness tem- 
perature 3D map Tskyia, 5, v) by applying the frequency or wave- 



Comparison P(k) LSS/21cm , Radio foregrounds 




1e.09 I 1 — ' — ' — I — 1-- 

0.02 0.04 0.06 0.1 0.15 0.2 0.4 

k_comov (h Mpc'^-I) 



Fig. 9. Comparison of the 21cm LSS power spectrum (red, 
orange) with the radio foreground power spectrum. The ra- 
dio sky power spectrum is shown for the GSM (Model-I) sky 
model (dark blue), as well as for our simple model based on 
Haslam-i-NVSS (Model-II, black). The curves with circle mark- 
ers show the power spectrum as observed by a perfect instrument 
with a 25 arcmin (FWHM) gaussian beam. 

length dependent instrument response *??(«, v, A). We have con- 
sidered the simple case where the instrument response is con- 
stant throughout the survey area, or independent of the sky di- 
rection. For each frequency Vk or wavelength - cjvk: 

1. Apply a 2D Fourier transform to compute sky angular 
Fourier amplitudes 

Tskyia, 6, Ak) ^ 2D - FFT T,ky{u, v, Ak) 

2. Apply instrument response in the angular wave mode 
plane. We use here the normalized instrument response 
'Riu,v,Ak) < 1. 

Tsky(u, V, Ak) — » Tskyiu, V, Ak) X Hiu, v, Ak) 

3. Apply inverse 2D Fourier transform to compute the mea- 
sured sky brightness temperature map, without instrumental 
(electronic/rjv.s) white noise: 

T,ky(u, V, Ak)x'R(u, V, A) Inv - 2D - FFT ^ T„,,,i(a, 6, Ak) 

4. Add white noise (gaussian fluctuations) to the pixel map 
temperatures to obtain the measured sky brightness temper- 
ature Tinesifl, 5, Vk). We have also considered that the system 
temperature and thus the additive white noise level was in- 
dependent of the frequency or wavelength. 

The LSS signal extraction depends indeed on the white noise 
level. The results shown here correspond to the (a) instrument 
configuration, a packed array of 11x11 = 121 dishes (5 meter 
diameter), with a white noise level corresponding to cr„„,5p = 
0.25 mK per 3 x 3arcmin^ x 500 kHz cell. 

A brief description of the simple component separation pro- 
cedure that we have applied is given here: 
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Haslam+NVSS@884 MHz*Lobe(35arcmin) 

Fig. 8. Comparison of GSM map (top) and Model-II sky map at 884 MHz (bottom). The Model-II (Haslam+NVSS) has been 
smoothed with a 35 arcmin (FWHM) gaussian beam. 



1. The measured sky brightness temperature is first corrected 
for the frequency dependent beam effects through a convolu- 
tion by a fiducial frequency independent beam. This correc- 
tion corresponds to a smearing or degradation of the angular 
resolution. We assume that we have a perfect knowledge of 
the intrinsic instrument response, up to a threshold numeri- 
cal level of about > 1% for 'R(u,v,A). We recall that this is 
the normalized instrument response, :R(m, v, A) < 1. 

T,„,Aa,6,v)^Ti{?;\a,5,v) 

The virtual target instrument has a beam width larger than 
the worst real instrument beam, i.e at the lowest observed 
frequency. 

2. For each sky direction (a, 5), a power law T - ^o^^) is 
fitted to the beam-corrected brightness temperature. The fit 
is done through a linear fit in the logio(r),logio(v) plane 
and we show here the results for a pure power law (PI) or 
modified power law (P2): 

PI : logio(r^:7(y)) = fl + /7logio(v/yo) 

P2 : logio(r^™») = fl + ^logi()(v/yo) + clogio(v/vo)2 

where b is the power law index and Tq = 10" is the brightness 
temperature at the reference frequency vq: 

3. The difference between the beam-corrected sky temperature 
and the fitted power law {Ti)(a, 6), b{a, 5)) is our extracted 21 
cm LSS signal. 



Figure 10 shows the performance of this procedure at a red- 



shift ~ 0.6, for the two radio sky models used here: GSM/Model- 
1 and Haslam-HNVSS/Model-11. The 21 cm LSS power spec- 
trum, as seen by a perfect instrument with a 25 arcmin (FWHM) 
gaussian frequency independent beam is shown in orange (solid 
line), and the extracted power spectrum, after beam correction 
and foreground separation with second order polynomial fit (P2) 
is shown in red (circle markers). We have also represented the 
obtained power spectrum without applying the beam correction 
(step 1 above), or with the first order polynomial fit (PI). 



Figure 11 shows a comparison of the original 21 cm bright- 
ness temperature map at 884 MHz with the recovered 21 cm 
map, after subtraction of the radio continuum component. It can 
be seen that structures present in the original map have been 
correctly recovered, although the amplitude of the temperature 
fluctuations on the recovered map is significantly smaller (factor 
~ 5) than in the original map. This is mostly due to the damping 
of the large scale {k < 0.04/!Mpc"') due the poor interferometer 
response at large angle {9 > 4°). 

We have shown that it should be possible to measure the red 
shifted 21 cm emission fluctuations in the presence of the strong 
radio continuum signal, provided that this latter has a smooth 
frequency dependence. However, a rather precise knowledge of 
the instrument beam and the beam correction or smearing pro- 
cedure described here are key ingredient for recovering the 21 
cm LSS power spectrum. It is also important to note that while 
it is enough to correct the beam to the lowest resolution instru- 
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ment beam (~ 30' or D ~ 50 meter @ 820 MHz) for the GSM 
sky model, a stronger beam correction has to be applied ((~ 36' 
or D ~ 40 meter @ 820 MHz) for the Model-II to reduce sig- 
nificantly the ripples from bright radio sources. We have also 
applied the same procedure to simulate observations and LSS 
signal extraction for an instrument with a frequency dependent 
gaussian beam shape. The mode mixing effect is greatly reduced 
for such a smooth beam, compared to the more complex instru 
ment response v. A) used for the results shown in figure 10 



TF(k): Extracted_P21(k)/P21(k) z=a.5,1.a,1.5,2.a,2.5 (blue...brawn) 



4.3. Piiik) measurement transfer function 

The recovered red shifted 21 cm emission power spectrum 
P^i(k) suffers a number of distortions, mostly damping, com- 
pared to the original Piiik) due to the instrument response and 
the component separation procedure. We expect damping at 
small scales, or larges k, due to the finite instrument size, but 
also at large scales, small k, if total power measurements (auto- 
correlations) are not used in the case of interferometers. The 
sky reconstruction and the component separation introduce addi- 
tional filtering and distortions. Ideally, one has to define a power 
spectrum measurement response or transfer function in the ra- 
dial direction, (A or redshift, T(A;||)) and in the transverse plane ( 
T(fcj_) ). The real transverse plane transfer function might even 
be anisotropic. 

However, in the scope of the present study, we define an 
overall transfer function T{k) as the ratio of the recovered 3D 
power spectrum Pj^X^) the original P2\{k): 



T{k) = p',r(k)/P2i(k) 



(35) 



Figure 12 shows this overall transfer function for the sim- 
ulations and component separation performed here, around 
z ~ 0.6, for the instrumental setup (a), a filled array of 121 
Ddish = 5 m dishes. The orange/yellow curve shows the ratio 
P^2i°'"'""'{k)/P2\(k) of the computed to the original power spec- 
trum, if the original LSS temperature cube is smoothed by the 
frequency independent target beam FWHM=30' for the GSM 
simulations (left), 36' for Model-II (right). This orange/yellow 
curve shows the damping effect due to the finite instrument 
size at small scales (k > 0.1 /i Mpc"', 6* < 1°). The recovered 
power spectrum suffers also significant damping at large scales 
k < 0.05 /iMpc"', due to poor interferometer response at large 
angles (6 > 4° -5°), as well as to the filtering of radial or longitu- 
dinal Fourier modes along the frequency or redshift direction (fey) 
by the component separation algorithm. The red curve shows 
the ratio of P{k) computed on the recovered or extracted 21 cm 
LSS signal, to the original LSS temperature cube P'^^{k)IP2i{k) 
and corresponds to the transfer function T{k) defined above, for 
z - 0.6 and instrument setup (a). The black (thin line) curve 
shows the ratio of recovered to the smoothed power spectrum 
P''^^'{k)IP'^l'""''"^(k). This latter ratio (black curve) exceeds one 
for k > 0.2, which is due to the noise or system temperature. It 
should be stressed that the simulations presented in this section 
were focused on the study of the radio foreground effects and 
have been carried intently with a very low instrumental noise 
level of 0.25 mK per pixel, corresponding to several years of 
continuous observations (~ 10 hours per 3' x 3' pixel). 

This transfer function is well represented by the analytical 
form: 

T«)=^xexp(-1) ,36) 

We have performed simulation of observations and radio 
foreground subtraction using the procedure described here for 




0.01 0.02 



0.1 0.2 0.5 

k_comov (h Mpc^-1) 



Fig. 13. Fitted/smoothed transfer function T(k) obtained for the 
recovered 21 cm power spectrum at different redshifts, z = 
0.5, 1.0, 1.5, 2.0, 2.5 for the instrument configuration (e), 20 x 20 
packed array interferometer. 



different redshifts and instrument configurations, in particular 
for the (e) configuration with 400 five-meter dishes. As the syn- 
chrotron and radio source strength increases quickly with de- 
creasing frequency, we have seen that recovering the 21 cm 
LSS signal becomes difficult for larger redshifts, in particular 
for z>2. 

We have determined the transfer function parameters of eq. 
[36]^A; kg, kc for setup (e) for three redshifts, z = 0.5, 1, 1.5, and 
then extrapolated the value of the parameters for redshift z = 
2, 2.5. The value of the parameters are grouped in table |5] and 
the smoothed transfer functions are shown on figure [13] 



z 


0.5 


1.0 


1.5 


2.0 


2.5 


kA 


0.006 


0.005 


0.004 


0.0035 


0.003 


ks 


0.038 


0.019 


0.012 


0.0093 


0.008 


kc 


0.16 


0.08 


0.05 


0.038 


0.032 



Table 5. Value of the parameters for the transfer function (eq. 
[36] l at different redshift for instrumental setup (e), 20x20 packed 
array interferometer. 



5. Sensitivity to cosmological parameters 

The impact of the various telescope configurations on the sen- 
sitivity for 21 cm power spectrum measurement has been dis- 
cussed in section |3] Fig. |6] shows the noise power spectra, and 
allows us to rank visually the configurations in terms of instru- 
ment noise contribution to P(k) measurement. The differences in 
Pnoise will translate into differing precisions in the reconstruction 
of the BAO peak positions and in the estimation of cosmological 
parameters. In addition, we have seen (sec. |4.2[ ) that subtraction 
of continuum radio emissions. Galactic synchrotron and radio 
sources, has also an effect on the measured 21 cm power spec- 
trum. In this paragraph, we present our method and the results 
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Fig. 10. Recovered power spectrum of the 21cm LSS temperature fluctuations, separated from the continuum radio emissions at 
z ~ 0.6, for the instrument configuration (a), 11x11 packed array interferometer Left: GSM/Model-I , right: Haslam+NVSS/Model- 
IL black curve shows the residual after foreground subtraction, corresponding to the 21 cm signal, WITHOUT applying the beam 
correction. Red curve shows the recovered 21 cm signal power spectrum, for P2 type fit of the frequency dependence of the radio 
continuum, and violet curve is the PI fit (see text). The orange/yellow curve shows the original 21 cm signal power spectrum, 
smoothed with a perfect, frequency independent gaussian beam. 



for the precisions on the estimation of Dark Energy parameters, 
through a radio survey of the redshifted 21 cm emission of LSS, 
with an instrumental setup similar to the (e) configuration (sec. 
|3.3| l, 400 five-meter diameter dishes, arranged into a filled 20x20 
array. 

5.1. BAO peak precision 

In order to estimate the precision with which BAO peak 
positions can be measured, we used a method similar to 
the one established in (Bla ke and Glazebro ok (2003)) and 
( |Glazebrookand Blake (2005)| ). 

To this end, we generated reconstructed power spectra 
P'''^'^(k) for slices of Universe with a quarter-sky coverage and 
a redshift depth, Az = 0.5 for 0.25 < z < 2.75. The peaks in the 
generated spectra were then determined by a fitting procedure 
and the reconstructed peak positions compared with the gener- 
ated peak positions. The reconstructed power spectrum used in 
the simulation is the sum of the expected Hi signal term, corre- 
sponding to equatio ns [9l and 10 damped by the transfer function 
T(k) (Eq. 36 , table ISjl and a white noise component P„oise cal- 
culated according to the equation 29 established in section 3.3 
with = 400: 



P'-"'\k)^P2i(k)xT{k) + P„ 



(37) 



where the different terms (P2\{k),T{k),P„oise) depend on the 
slice redshift. The expected 21 cm power spectrum P2\{k) has 



been generated according to the formula: 

— - \ + 

Pref(kj_,k\{) 

f ~ 

Akexp{-(k/T)")sm 2n ^ - 
\ k 



k^ 

ii_ 

2 k^ 
BAO± '^BAOW ) 



(38) 



where k - ^k\ + the parameters A, a and r are adjusted to 

the formula presented in ( |Eisenstein & Hu (1998)| l. Pre]iki_,k\\} 
is the envelop curve of the HI power spectrum without baryonic 
oscillations. The parameters A;bao± and kgj^ow are the inverses 
of the oscillation periods in k-space. The following values have 
been used for these parameters for the results presented here: 
A - 1.0, T = 0.1 /iMpc"', a - 1.4 and ^bao± = ^bao|| = 
0.060 /zMpc-i. 

Each simulation is performed for a given set of param- 
eters which are: the system temperature,?",,.,, an observation 
time, tohs, an average redshift and a redshift depth, Az = 0.5. 
Then, each simulated power spectrum is fitted with a two di- 
mensional normalized function Ptoi(k±, k\\)/Pref{k±,k\{) which is 
the sum of the signal power spectrum damped by the trans- 
fer function and the noise power spectrum multiplied by a lin- 
ear term, ao + aik. The upper limit k,nax in k of the fit cor- 
responds to the approximate position of the Unear/non-linear 
transition. This limit is estabUshed on the basis of the crite- 
rion discussed in ( |Blake and Glazebrook (2003)1 ). In practice. 
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Fig. 11. Comparison of the original 21 cm LSS temperature map @ 884 MHz (z ~ 0.6), smoothed with 25 arcmin (FWHM) beam 
(top), and the recovered LSS map, after foreground subtraction for Model-I (GSM) (bottom), for the instrument configuration (a), 
11x11 packed array interferometer. Notice the difference between the temperature color scales (mK) for the top and bottom maps. 



we used for the redshifts 



0.5, 1.0 and 1.5 respectively 



0.145/iMpc" 



0.18 /iMpc"' and 0.23 /jMpc"'. 



Figure 14 shows the result of the fit for one of these simu- 
lations. Figure[T5]histograms the recovered values of kBAOi. and 
ksAOW f^i" 100 simulations. The widths of the two distributions 
give an estimate of the statistical errors. 

In addition, in the fitting procedure, both the parameters 
modeling the signal A, r, a and the parameter correcting the 
noise power spectrum (flo,ai) are floated to take into account 
the possible ignorance of the signal shape and the uncertainties 
in the computation of the noise power spectrum. In this way, we 
can correct possible imperfections and the systematic uncertain- 
ties are directly propagated to statistical errors on the relevant 
parameters ksAOi. and ksAow - By subtracting the fitted noise con- 
tribution to each simulation, the baryonic oscillations are clearly 



observed, for instance, on Fig. 16 



In our comparison of the various configurations, we have 
considered the following cases for Az = 0.5 slices with 0.25 < 
z < 2.75. 



- Simulation without electronics noise: the statistical errors on 
the power spectrum are directly related to the number of 
modes in the surveyed volume V corresponding to Az = 0.5 
slice with the solid angle Q.,ot = 1 tt sr. The number of mode 



z 


0.5 


1.0 


1.5 


2.0 


2.5 


P„o,sc mK- (Mpc/h)^ 


8.5 


35 


75 


120 


170 



Table 6. Instrument or electronic noise spectral power P„oise for 
a N - 400 dish interferometer with T^ys = 50 K and = 1 
year to survey f2,o/ = r: si 



Nsk in the wave number interval 6k can be written as: 



V 



k^6k 



(39) 



Noise: we add the instrument noise as a constant term P„ 



as described in Eq. 29 Table |6] gives the white noise level 
for Tsys - 50K and one year total observation time to survey 
Q,„, = 1 TT sr. 

- Noise with transfer function: we take into account of the in- 
terferometer and radio foreground subtraction represented as 
the measured P(k) transfer function T{k) (section |4.3| l, as 
well as instrument noise P„oise- 

Table|7]summarizes the result. The errors both on kBAOi. and 
ksAOW decrease as a function of redshift for simulations without 
electronic noise because the volume of the universe probed is 
larger Once we apply the electronics noise, each slice in red- 
shift give comparable results. Finally, after applying the full re- 
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Fig. 12. Ratio of the reconstructed or extracted 2 1 cm power spectrum, after foreground removal, to the initial 21cm power spectrum, 
T(k) = P2[(k)/P2i(k), at z ~ 0.6, for the instrument configuration (a), 11x11 packed array interferometer. Left: GSM/Model-I , 
right: Haslam+NVSS/Model-II. 
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k (h Mpc ') 

Fig. 14. ID projection of the power spectrum for one simulation. 
The Hi power spectrum is divided by an envelop curve P{k)ref 
corresponding to the power spectrum without baryonic oscilla- 
tions. The dots represents one simulation for a "packed" array of 
cylinders with a system temperature, Tj,,, = 50K, an observation 
time, Tohs - 1 year, a solid angle of Ijrsr, an average redshift, 
z = 1.5 and a redshift depth, Az = 0.5. The solid line is the result 
of the fit to the data. 



construction of the interferometer, the best accuracy is obtained 
for the first slices in redshift around 0.5 and 1 .0 for an identical 
time of observation. We can optimize the survey by using a dif- 
ferent observation time for each slice in redshift. Finally, for a 3 
year survey we can split in five observation periods with dura- 
tions which are 3 months, 3 months, 6 months, 1 year and 1 year 
respectively for redshift 0.5, 1.0, 1.5, 2.0 and 2.5. 



Mean 0.06038 
RMS 0.003212 




Fig. 15. Distributions of the reconstructed wavelength ksAOi. 
and ksAOW respectively, perpendicular and parallel to the line of 
sight for simulations as in Fig. 14 The fit by a Gaussian of the 
distribution (solid line) gives the width of the distribution which 
represents the statistical error expected on these parameters. 



5.2. Expected sensitivity on wq and Wg 

The observations give the Hj power spectrum in angle-angle- 
redshift space rather than in real space. The inverse of the peak 
positions in the observed power spectrum therefore gives the an- 
gular and redshift intervals corresponding to the sonic horizon. 
The peaks in the angular spectrum are proportional to dT{z)/as 
and those in the redshift spectrum to dH(z)/as. ~ 105/z^'Mpc 
is the acoustic horizon comoving size at recombination, driz) = 
(1 + z)dA is the comoving angular distance and dH - c/H{z) is 
the Hubble distance (see Eq.|6]l: 



dH = 



c/Ho 



dr ^ I dn 
Jo 



{z)dz (40) 
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2.8 


(3-months/redshift) 


Cr{kBAO\\)lkBAO\\ {%) 


4.1 
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4.3 
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Noise witli Transfer Function 


'r(kBAO±)/kBAO± (%) 


3.0 


2.5 


3.5 


5.2 


6.5 


(3-months/redshift) 


a-(kBAOll)/kBAO\l (%) 


4.8 


4.0 


6.2 


9.3 


10.3 


Optimized survey 


Cr(kBAO±)/kBAO± (%) 


3.0 


2.5 


2.3 


2.0 


2.7 


(Observation time : 3 years) 


(r{kBAOll)/kBAO\l (%) 


4.8 


4.0 


4.1 


3.6 


4.3 



Table 7. Sensitivity on the measurement of kg^o^ ™d ^bao|| a function of the redshift z for various simulation configuration. 
row: simulations without noise with pure cosmic variance; 2"'' row: simulations with electronics noise for a telescope with dishes; 
3* row: simulations with same electronics noise and with correction with the transfer function ; 4* row: optimized survey with a 
total observation time of 3 years (3 months, 3 months, 6 months, 1 year and 1 year respectively for redshift 0.5, 1.0, 1.5,2.0 and 2.5 
). 
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Fig. 16. ID projection of the power spectrum averaged over 100 
simulations of the packed cylinder array b. The simulations are 
performed for the following conditions: a system temperature, 
Tsys = 50K, an observation time, T„i,s = 1 year, a solid angle of 
Insr, an average redshift, z = 1.5 and a redshift depth, Az - 0.5. 
The Hi power spectrum is divided by an envelop curve P{k)ref 
corresponding to the power spectrum without baryonic oscilla- 
tions and the background estimated by a fit is subtracted. The 
errors are the RMS of the 100 distributions for each k bin and 
the dots are the mean of the distribution for each k bin. 
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Fig. 17. The two "Hubble diagrams" for BAO experiments. The 
four falling curves give the angular size of the acoustic hori- 
zon (left scale) and the four rising curves give the redshift inter- 
val of the acoustic horizon (right scale). The solid lines are for 
{Qm,^a, w) = (0.27,0.73,-1), the dashed for (1,0,-1) the dot- 
ted for (0.27,0,-1), and the dash-dotted for (0.27,0.73,-0.9), 
The error bars on the solid curve correspond to the four-month 
run (packed array) of Table |7] 



The quantities dj, dn and a, all depend on the cosmological pa- 
rameters. Figure 17 gives the angular and redshift intervals as 



a function of redshift for four cosmological models. The eiTor 
bars on the lines for (Qjv/.^^a) = (0.27,0.73) correspond to the 
expected errors on the peak positions taken from Table|7]for the 
four-month runs with the packed array. We see that with these 
uncertainties, the data would be able to measure w at better than 
the 10% level. 

To estimate the sensitivity to parameters describing dark en- 
ergy equation of state, we follow the procedure explained in 
( |Blake and Glazebrook (2003)] ). We can introduce the equation 
of state of dark energy, w(z) - wq + Wa ■ zl (_l + z) replacing 
□a in the definition of dj(z) and djj{z), (Eq. 40 1 by: 



Qa — > exp 



, r 1 + w(z') 

Jo l+z' 



dz 



(41) 



where is the present-day dark energy fraction with respect to 
the critical density. Using the relative errors on kgj^o^ ™d k^Aow 
given in Tab. [7] we can compute the Fisher matrix for five cos- 
mological parameter: (Q,„, Q;,, h, wo, Wa)- Then, the combination 



of this BAO Fisher matrix with the Fisher matrix obtained for 
Planck mission, allows us to compute the errors on dark energy 
parameters. The Planck Fisher matrix is obtained for the 8 pa- 
rameters (assuming a flat universe): Q,„, Q;,, h, wq, Wg, erg, 
(spectral index of the primordial power spectrum) and r (optical 
depth to the last-scatter surface). 

For an optimized project over a redshift range, 0.25 < z < 
2.75, with a total observation time of 3 years, the packed 400- 
dish interferometer array has a precision of 12% on wq and 48% 
on Wa- The Figure of Merit, the inverse of the area in the 95% 
confidence level contours is 38. Finally, Fig. 18 shows a compar- 
ison of different BAO projects, with a set of priors on (f2„,, Q^,, h) 
coiTesponding to the expected precision on these parameters in 
early 2010's. This BAO project based on Hj intensity mapping 
is clearly competitive with the current generation of optical sur- 
veys such as SDSS-III ( 1SDSS-III(20()8)] |. 
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Fig. 18. lo" and 2cr confidence level contours in the parameter 
plane (wq, Wa) for two BAO projects: SDSS-Ill (LRG) project 
(blue dotted line), 21 cm project with HI intensity mapping 
(black solid line). 



6. Conclusions 

The 3D mapping of redshifted 21 cm emission though Intensity 
Mapping is a novel and complementary approach to optical sur- 
veys to study the statistical properties of the large scale struc- 
tures in the universe up to redshifts z < 3. A radio instrument 
with large instantaneous field of view (10-100 deg^) and large 
bandwidth (> 100 MHz) with ~ 10 arcmin resolution is needed 
to perform a cosmological neutral hydrogen survey over a sig- 
nificant fraction of the sky. We have shown that a nearly packed 
interferometer array with few hundred receiver elements spread 
over an hectare or a hundred beam focal plane array with a 
~ 100 meter primary reflector will have the required sensitiv- 
ity to measure the 21 cm power spectrum. A method to com- 
pute the instrument response for interferometers has been de- 
veloped and we have computed the noise power spectrum for 
various telescope configurations. The Galactic synchrotron and 
radio sources are a thousand time brighter than the redshifted 21 
cm signal, making the measurement of this latter signal a ma- 
jor scientific and technical challenge. We have also studied the 
performance of a simple foreground subtraction method through 
realistic models of the sky emissions in the GHz domain and 
simulation of interferometric observations. We have been able 
to show that the cosmological 21 cm signal from the LSS should 
be observable, but requkes a very good knowledge of the instru- 
ment response. Our method has allowed us to define and com- 
pute the overall transfer function or response function for the 
measurement of the 21 cm power spectrum. Finally, we have 
used the computed noise power spectrum and P{k) measurement 
response function to estimate the precision on the determination 
of Dark Energy parameters, for a 21 cm BAO survey. Such a 
radio survey could be carried using the current technology and 
would be competitive with the ongoing or planned optical sur- 
veys for dark energy, with a fraction of their cost. 
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